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Abstract 

Using the the convex semidefinite programming method and superoperator formal- 
ism we obtain the finite quantum tomography of some mixed quantum states such 
as: qudit tomography, N-qubit tomography, phase tomography and coherent spin state 
tomography, where that obtained results are in agreement with those of References 
[21, 24, 25, 4, 26]. Keywords: finite quantum tomograpy, Semi-definite program- 
ming,superoperator formalism, qubit quantum tomography and truncation. 
PACs Index: 03.65.Ud 
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The quantum complementarity principle does not allow to recover the quantum state from 
measurements on a single system, unless we have some prior information on it. On the other 
hand, the no cloning theorem ensures that it is not possible to make exact copies of a quantum 
system, without having prior knowledge of its state. Hence, the only possibility for devising 
a state reconstruction procedure is to provide a measuring strategy that employs numerous 
identical (although unknown) copies of the system, so that different measurements may be 
performed on each of the copies. 

The problem of state estimation resorts essentially to estimating arbitrary operators of a 
quantum system by using the result of measurements of a set of observables. If this set of 
observables is sufficient to give full knowledge of the system state, then we define it a quorum. 
Notice that, in general, a system may allow various, different quorums. Quantum tomography 
was born [1, 2] as a state reconstruction technique in the optical domain, and has recently been 
extended [3] to a vast class of systems. By extension, we now denote as Quantum Tomography 
all unbiased quantum state reconstruction procedures, i.e. those procedures which are affected 
only by statistical errors that can be made arbitrarily small by increasing the number of 
measurements. Tomography makes use of the results of the quorum measurements in order to 
reconstruct the expectation value of arbitrary operators (even not observables) acting on the 
system Hilbert space. 

In principle, a precise knowledge of the density matrix would require an infinite number 
of measurements on identical preparations of radiation. However, in real experiments one 
has only a finite number of data at ones disposal, and thus a statistical analysis and errors 
estimation are needed. 

Authors of Ref. [4] presented several schemes for a reconstruction of states of quantum 
systems from measured data: 
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(1) The maximum entropy (MaxEnt) principle leads to a complete reconstruction of quan- 
tum states, i.e. quantum states are uniquely determined. 

(2) Quantum systems can be estimated with the help of quantum Bayesian inference. 

(3) Estimation of a quantum state with the highest fidelity and Showed how this optimal 
measurement can in principle be realized [4]. 

On the other hand, over the past years, semidefinite programming (SDP) has been rec- 
ognized as valuable numerical tools for control system analysis and design. In (SDP) one 
minimizes a linear function subject to the constraint that an affine combination of symmetric 
matrices is positive semidefinite. SDP, has been studied (under various names) as far back as 
the 1940s. Subsequent research in semidefinite programming during the 1990s was driven by 
applications in combinatorial optimization [5], communications and signal processing [6, 7, 8], 
and other areas of engineering [9]. Although semidefinite programming is designed to be ap- 
plied in numerical methods it can be used for analytic computations, too. Some authors try to 
use the SDP to construct an explicit entanglement witness [10, 11]. Kitaev used semidefinite 
programming duality to prove the impossibility of quantum coin flipping [12], and Rains gave 
bounds on distillable entanglement using semidefinite programming [13]. In the context of 
quantum computation, Barnum, Saks and Szegedy reformulated quantum query complexity 
in terms of a semidefinite program [14]. The problem of finding the optimal measurement to 
distinguish between a set of quantum states was first formulated as a semidefinite program 
in 1972 by Holevo, who gave optimality conditions equivalent to the complementary slackness 
conditions [15]. Recently, Eldar, Megretski and Verghese showed that the optimal measure- 
ments can be found efficiently by solving the dual followed by the use of linear programming 
[16]. Also in [17] used semidefinite programming to show that the standard algorithm imple- 
ments the optimal set of measurements. All of the above mentioned applications indicate that 
the method of SDP is very useful. 

In a laboratory and in practice, we always deal with finite ensembles of copies of the 
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measured system. This implies the need of developing novel tools specially designed to process 
realistic and finite experimental samples. Then it is necessary to truncate the Hilbert space 
to a finite dimensional basis [18]. In this paper we use the SDP method in order to obtain 
quantum tomography with truncating the infinite Banach space to a finite dimensional basis. 

The paper is organized as follows: 
In section-2 we define semidefinite programming. In section -3 we define superoperator formal- 
ism. In section -4 we describe the projection method and using SDP method and superoperator 
formalism we obtain finite quantum tomography. In section -5 we obtain some typical finite 
quantum tomographic examples, such as: finite dimensional qudit quantum tomography N- 
qubit tomography finite dimensional phase tomography and coherent spin state tomography 
with SDP method and superoperator formalism. The paper is ended with a brief conclusion. 

2 Semi-definite programming 

A SDP is a particular type of convex optimization problem [19]. A SDP problem requires 
minimizing a linear function subject to a linear matrix inequality (LMI) constraint [20]: 

minimize V = c T x 

(2-1) 

subject to F(x) > 0, 

where c is a given vector, x T = (xi, ...,x n ), and F(x) = F + J2i x iFi, for some fixed hermitian 
matrices Fj. The inequality sign in F(x) > means that F(x) is positive semidefinite. 

This problem is called the primal problem. Vectors x whose components are the variables 
of the problem and satisfy the constraint F(x) > are called primal feasible points, and if 
they satisfy F(x) > they are called strictly feasible points. The minimal objective value c T x 
is by convention denoted as V* and is called the primal optimal value. 

Due to the convexity of set of feasible points, SDP has a nice duality structure, with, the 
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associated dual program being: 

maximize —Tr[F Z] 

Z > (2-2) 
Tr[F t Z] = a. 

Here the variable is the real symmetric (or Hermitean) matrix Z, and the data c, Fj are 
the same as in the primal problem. Correspondingly, matrices Z satisfying the constraints are 
called dual feasible (or strictly dual feasible if Z > 0). The maximal objective value —TrF Z, 
the dual optimal value, is denoted as d*. 

The objective value of a primal(dual) feasible point is an upper (lower) bound on V*(d*. 
The main reason why one is interested in the dual problem is that one can prove that d* <V*, 
and under relatively mild assumptions, we can have V* = d*. if the equality holds, one can 
prove the following optimality condition on x: 

A primal feasible x and a dual feasible Z are optimal which is denoted by x and Z if and 
only if 

F(x)Z = ZF(x) = 0. (2-3) 

This latter condition is called the complementary slackness condition. 

In one way or another, numerical methods for solving SDP problems always exploit the 
inequality d < d* < V* < V, where d and V are the objective values for any dual feasible 
point and primal feasible point, respectively. The difference 

V - d = c T x + Tr[F Z] = Tr[F(x)Z] > (2-4) 

is called the duality gap. If the equality d* = V* holds, i.e., the optimal duality gap is zero, 
then we say that strong duality holds. 
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3 Superoperator formalism 

In order to treat discrete and continuous density operator representations on an equal footing, 
we introduce the following superoperator formalism. The set of linear operators acting on a D- 
dimensional Hilbert space H is a D 2 -dimensional complex vector space C(H). Let us introduce 
operator "kets" | A) — A and "bras" (A \— A', distinguished from vector kets and bras by the 
use of round brackets. Then the natural inner product on L(H), the trace-norm inner product, 
can be written as (A | B) — tr(A^B). The notation S —\ A)(B \ defines a superoperator S 
acting like 

S\X)=\ A)(B | X) = tr(B*X)A. (3-5) 

Now let the set {| Nj)} constitute a (complete or overcomplete) operator basis; i.e., let the 
operator kets | Nj) span the vector space L(H). It follows that the superoperator Q defined by 

= E I N j)( N i I ( 3 - 6 ) 

3 

is invertible. The operators 

Qj = I Nj) (3-7) 
form a dual basis, which gives rise to the following resolutions of the superoperator identity: 

i = £l^)(^h£l^)(Q,-|. (3-8) 
j j 

An arbitrary operator A can be expanded as 

A = E I Ni)(Qi I A )=H NMQjA) ( 3 - 9 ) 
j j 

and 

^ = E I QiW | A) = Y.QMN}A) (3-10) 

j j 

These expansions are unique if and only if the operators Nj are linearly independent [21]. 
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4 Projection method as a semidefinite programming and 
finite quantum tomography 



4.1 Bases and frames 

In this section we collect some rudimentary facts that will be used in what follows. 
A basis is one of the most fundamental concepts in linear algebra. 

A set of linearly independent vectors {ej}™ =1 in a finite dimensional complex vector space 
V is a basis for V if, for each / e V, there exist coefficients ci, C2, c n G C such that 

n 

f = J2 c i e i- ( 4 - n ) 

1=1 

The independence condition implies that the coefficients c±, ...,c n are unique. 

For infinite dimensional vector spaces, the concept of a basis is more complicated. 
An {ejjg^ C H is an orthonormal system (ONS) [22] if 

<e i ,e j >=5 ij . (4-12) 
An ONS {e i }'^ l is an orthonormal basis (ONB) if 

U = span{ei\Zi (4-13) 
when {ei}^ is an ONB, each / e H can be written as 

oo 

f = J2<f,e l >e t . (4-14) 
i=i 

Definition: Two sequences {a^} and {yj} in a Hilbert space Ti are said to be biorthonormal, 

if 

<x i ,y j >=5 ij . (4-15) 

A sequence {y^} biorthogonal to a basis {xj} for 7i is itself a basis for H, and we have for 
each x the representation 

oo oo 

x = ^2 < x,Ui > Xi, and a; = ^ < x, x» > y*. (4-16) 

i=l i=l 
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Frame: A family of elements {fi}iei Q H is called a frame for H if there exist constants 
A, B > such that 

A\\f\\ 2 < £ | < fji > | 2 < 5||/|| 2 , V/ e H, (4-17) 

where / is a countable index set. The numbers A, B are called frame bounds. They are 
not unique. The optimal frame bounds are the biggest possible value for A and the smallest 
possible value for B in (4-17). If we can choose A = B, the frame is called tight. If a frame 
ceases to be a frame when any element is removed, the frame is said to be exact. Since a frame 
{fi}iei is a Bessel sequence, the operator 

T : Z 2 (J) - H , T{ Cl ieI } = £ ah (4-18) 

iei 

is bounded and linear; T is sometimes called the preframe operator. The adjoint operator is 
given by 

T* : H - l 2 (I) , T*f = {< /, ft >}Zv (4-19) 
By composing the operators T and T*, we obtain the operator 

oo 

S-.n^H, Sf = TT*f = ]T< f,fi> & (4-20) 

1=1 

where S is called the frame operator with 

AI < S < BI. (4-21) 
The frame operator is a bounded, positive, and invertible operator. 

4.2 Frames in finite-dimensional spaces 

We investigate the properties of a frame generated by a finite subset of a Hilbert space. 

Calculation of the frame coefficients {< /, S* -1 ^ >} involves inversion of the frame operator 
S. In practice it can be a problem if the underlying Hilbert space is infinite dimensional. There 
is an approach to the problem as follows [22] : 
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Given the frame we consider finite subsets {fi}2 = i,n G N. It can be shown that 

{fi}i=i is a frame for H n = span{/j}™ and the corresponding frame operator is S n : H n — > 7i n 
and the orthogonal projection P n on 7i n is 

n 

Pnf = J2<f,S^ 1 f i >f i , feH. (4-22) 
1=1 

For n — > oo, P n / — > / = X^i < /, > /i, one can hope that the coefficients < 

fjS^fi >converges to the frame coefficients for /, i.e., that 

< /, S- 1 ^ >^< f, S-'f- > as n oo, e I, V/ G H. (4-23) 

If (4-23) is satisfied we say that the projection method works. In this case the frame coefficients 
can be approximated as close as we want using finite dimensional methods, i.e., linear algebra, 
since S n is an operator on the finite dimensional space Ti n . This is a very important property 
for applications: for example, it makes it possible to use computers to approximate the frame 
coefficients. 

In the following subsection we obtain finite quantum tomography using semidefinite pro- 
gramming. 

4.3 Finite quantum tomography via semidefinite programming 

Quantum state reconstruction schemes can be understood as an a posterior estimation of 
density operator of a given quantum mechanical system based on data obtained with the help 
of a macroscopic measurement apparatus. Only if an infinite ensemble is given can one find 
out the state. But infinite ensembles don't exist in practice. In a laboratory and in practice, 
we always deal with finite ensembles of copies of the measured system. This implies the need of 
developing novel tools specially designed to process realistic and finite experimental samples. 
Then it is necessary to truncate the Hilbert space to a finite dimensional basis [18]. 

Now in this work using the projection method and semidefinite programming we express the 
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mathematical structure correspond to finite tomography and obtain the tomographic formula 
based on finite Banach space. 

At first from (3-9)or (3-10) we assume that 

p = E I QAM I p) = E QMrfp) ( 4 - 24 ) 

j j 

is a density matrix in infinite dimensional Banach space, where {Nj} constitute a operator 
basis in superoperator formalism. Also let 

n 

P n = HW N i) (4-25) 

be a density matrix in finite dimensional banach space which is obtained from truncating the 
infinite dimensional Banach space. 

Using the properties of density matrix we have 

P-P n >0, (4-26) 

which in comparison with semidefinite programming we get 

F = p, Fj =| Nj) and Xj = Xj, for j = 1, ...,n. 

If we use the complementary slackness condition, and for a feasible (Z, A, m ), for j = 
1, n, we have 

Z(p - p n ) = 0, (4-27) 

or 

Z(p-jr\ 1 \N J ))=0. (4-28) 

3=1 

Using resolution of the superoperator identity (3-8) we obtain 

J2Z\N i )(Q t \[p-J2^\N j )] = 0, for j = l,..., n (4-29) 

» 3 

Therefore, we have 

]T(Z | NJKQ, | p) - \] = , i = 1, ...n. (4-30) 
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It is obvious that (Z \ Ni) = for i > n then we conclude that | iVj) G kerZ. Then we obtain 

*i = (Qi\p)=tr\pN}]. (4-31) 

Therefore we obtain the tomography formula in finite dimensional Banach space as the 
follow: 

n n 

P n = E I N i)(Qi I P) = E I Ni)tr[pN}}. (4-32) 
i=i i=i 

In the following, we will consider density matrix with orthogonal states of the form: 

oo 

p = EMpI^)(^I)I^X^I 

3 

where is a density matrix in infinite dimensional Hilbert space. In the superoperator formalism 
we can write 

\N i )=\Q j ) = \1> i )(il>i\ ( 4 - 33 ) 

Also let 

n n 

p n = E\-l^) = E^I^)(^l 

be a density matrix in finite dimensional Hilbert space which is obtained from truncating the 
infinite dimensional Hilbert space and is an orthogonal state. 

Using (4-32) we obtain the tomography formula in finite dimensional Hilbert space as the 
follow: 

n 

p n = EMp|^X^I)l^)(^l- (4-34) 
i=i 

In the following we describe some examples for finite dimensional quantum tomogarphy. 
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5 Some examples for finite quantum tomography with 
semidefinite programming 



5.1 Qudit tomography 

We begin with the set of Hermitian generators of SU(D); the generators, denoted by Xj , 
are labeled by a Roman index taken from the middle of the alphabet, which takes on values 
j = 1, ...,D 2 — 1 [23]. We represent the generators in an orthonormal basis \a >, labeled by 
a Roman letter taken from the beginning of the alphabet, which takes on values a = 1, D. 
With these conventions the generators are given by 

./ 1 I* 1: 

A, = T a = 1 (j: \b)(b\ - (a - l)\a)(a\) , 2 < a < D , (5-35) 
^a(a — 1) \6=i / 

j = D,...,(D + 2)(D-l)/2: 

A, = rit ) = i=(|a)(6| + |6)(a|) , 1 < a < b < D , (5-36) 
j = D(D + 2)/2,...,D 2 -l: 

A J = ri- ) = ^i(|a)(&|-|6)(a|) , 1 < a < b < D . (5-37) 

In Eqs.(5-36) and (5-37), the Roman index j stands for the pair of Roman indices, ab, whereas 
in Eq.(5-35), it stands for a single Roman index a. The generators are traceless and satisfy 

AjA fe = -jjbjk + djkAi + ifjkAi ■ (5-38) 

Here and wherever it is convenient throughout this paper, we use the summation convention 
to indicate a sum on repeated indices. The coefficients fjti, the structure constants of the 
Lie group SU(-D), are given by the commutators of the generators and are completely anti- 
symmetric in the three indices. The coefficients djki are given by the anti-commutators of the 
generators and are completely symmetric. 
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By supplementing the D 2 — 1 generators with the operator 

A = -±=I , (5-39) 
V D 



where I is the unit operator, we obtain a Hermitian operator basis for the space of linear 
operators in the qudit Hilbert space. This is an orthonormal basis, satisfying 

tr(A Q A /3 ) = 5 a(3 . (5-40) 

Here the Greek indices take on the values 0, . . . , D 2 — 1; throughout this paper, Greek indices 
take on D 2 or more values. Using this orthonormality relation, we can invert Eqs.(5-35)-(5-37) 
to give 



D 



\ a ^ a \ = 7i+ T= = [-(a-l)r,+ £ r 6 | , (5-41) 

D .ln( n _ u \ b=a+1 



\Ja(a- 1) 



l«><6| = ^(rff+ir^), 1 <a<b<D, (5-42) 

|6><«l = ^(r^-ir^), l<a<b<D. (5-43) 

Any qudit density operator can be expanded uniquely as 

1 

P = pC a \ a , (5-44) 
where the (real) expansion coefficients are given by 

c a = Dtr(p\ a ) . (5-45) 



Normalization implies that c = yD, so the density operator takes the form 

p = l(/ + Cj A,) = -i(/ + c-A). (5-46) 
Here c— Cjej can be regarded as a vector in a (D 2 — l)-dimensional real vector space, spanned 

— * 

by the orthonormal basis e}, and A = XjCj is an operator- valued vector. 
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In order to treat discrete density operator representation for a qudit we introduce the 
superoperator formalism and SDP method. Consider a discrete set of projection operators [23] 
define in finite dimensional Banach space 

N^=\n^><rf a \=^{l + ^.rQ , a = l,...,K. (5-47) 

The corresponding superoperator, 



where, orthonormal eigenoperators of Q are A = I/yD and T = J2j I \ >< \ I- 

We are now prepared to write the inverse of Q with respect to the left-right action as 

K \D + 1 D J K ' 

Thus the dual operators are given by 

\Q na ) = Q-'IKJ = D(D K +1 \ \ N a ) - ii-). (5-50) 

Using SDP method we get 

1 

F = — (1 + c.A) , f a =\N a ) and x a = A a for a = l,...,K. (5-51) 

From cmplementary slackness condition we have 

K = {Qn a \p). (5-52) 

Therefore, tomography relation in finite dimensional Banach space can be represented in the 
form 

K K D(D + 1) K 

P K = H I N na )(Q na | p) = £ Tr[Ql aP ]N na = 1 J 1 £ N a (I + Tr[N na p]). (5-53) 

a=l a=l a=l 

A qubit is two-level system, for which D = 2. There is a one-to-one correspondence between 
the pure states of a qubit and the points on the unit sphere, or Bloch sphere [21]. Any pure 
state of a qubit can be written in terms of the Pauli matrices (a±, <t 2 , <t 3 ), as 
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where rf = (ni; 712; 713) is a unit vector, and 1 denotes the unit matrix. An arbitrary state p, 
mixed or pure, of a qubit can be expressed as 

P=\(l + t.T?) (5-54) 

where <| S |< 1. 

In order to treat discrete density operator representation for a qubit we introduce the super- 
operator formalism and SDP method. Consider a discrete set of projection operators[21] in 
superoperator formalism 

=| n£ >< |= ^(1 + ~rf.h~t) , a = l,...,K. (5-55) 
The corresponding superoperator, 

= E I N n>J( N nZ l= -A K l !)(! I +Et^- 1 ^X 1 I + 1 I -^1+E l l EWMl. 

a=l a a 

(5-56) 

generates dual-basis operators and expansion coefficients proportional to those for the contin- 
uous representation [21] if and only if 

o = E^« ( 5 - 57 ) 

a 

When these conditions are satisfied, the superoperator (5-56) simplifies to 



3 , 



with an inverse 

G- 1 = 

K 



G- 1 = hm\+zY,\°^M (5-59) 



J 



which generates dual-basis operators 



Finite quantum tomography via semidefinite programming 17 

Then the density matrix in finite dimensional Banach space is given by (4-25). Using SDP 
method we get 

F = ^(1 + S.a) , f a =\N a ) and x a = \ a for a = l,...,K. (5-61) 

From cmplementary slackness condition we have 

\* = (Qn a | p) = Tr i^( l + 3 o>. r? Q )(l + S^) = -^(1 + 3~£.Jt a ))}. (5-62) 

Therefore, tomography relation (4-32) in finite dimensional Banach space can be represented 
in the form 

P K = E I Q a )(N a | p) = 1 E JV Q (1 + 3S^ Q ), (5-63) 

a=l ^ a=l 

For M qubits, we define the pure-product-state projector 

N(a) = N ai <g> ... <g> N aM = ^(1 + sn ai ) <g> ... <g> (1 + s.ra a J, (5-64) 

and 

Q(a)=Q nai <S)...<S)Q naM =— (l + 3s-ra ai )<g>...<g> (l + 3s-ra a J, (5-65) 

where n stands for the collection of unit vectors ni, ...,riM- Any M-qubit density operator can 
be expanded as 

P^ = E I W«)(Qa I P) = -L J2 N ai (l + 3^.lf ai )^...N aM (l + 3^.lt aM ), (5-66) 

where thus obtained result is in agreement with those of already obtained by one of the authors 
in [21, 4]. 



5.2 Phase tomography 

one possible means of describing the phase of a quantum mechanical fields is in terms of the 
Pegg-Barnett hermitian phase operator <3> [24, 25, 4]. This operator is defined in a finite (but 
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arbitrary large) dimensional Hilbert space. In a (s+l)-dimensional Hilbert space the phase 
state are defined as 



e >-- 



1 



J- m — n 



Vs + 1 n=0 

this Hilbert space is spanned by a complete orthonormal set of basis phase state | 6 m >, given 
by (5-67) with 

0-77" 777, 

O m = + —- , m = 0,l,...,s, (5-68) 
s + 1 

where #o is a reference phase. In terms of the state | 8 m > the Hermitian phase operator is 

s 

®e = E Q m I #m ><O m \. (5-69) 

m=0 

From the definition of the phase state (5-67), we can express the projector 9 m \ 9 m >< 9 m \ in 
terms of the number state basis: 

I 9 m >< 9 m |= (s + l)- 1 E E e*^'"")* | n' >< n | . (5-70) 

n=0 n'=0 

In this case $e is orthonormal then we can write the tomography using semidefinite pro- 
gramming. 

At first we assume that 

p = [ Tr(p$ e )®edLi8, (5-71) 
Je 

is a density matrix in infinite dimensional Banach space. Also let 

p> = J2^\®e), (5-72) 
e 

be a density matrix in finite dimensional Banach space which is obtained from truncating the 
infinite dimensional Banach space. Using the properties of density matrix we have 

P ~ p' > 0, (5-73) 

which is comparison with semidefinite programming we get 

F = p , F e =| $0) and x e = \e , for 6 = , 6 + 2tt. (5-74) 
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If we use the complementary slackness condition, and for a feasible (Z,\g max ), for 9 = 
9 , ...,9 + 2n, we have 

Z(p - p') = or Z(p - \ e | = 0. (5-75) 
Similar to superoperator formalism we obtain 

\ e = ($* | p) = Tr[p$ e }. (5-76) 
Therefore we obtain the tomography formula in finite dimensional Hilbert space as the follow: 

p' = El I P) = E I ^)Tr[p$ e ]. (5-77) 

If we generalized it when 9 is continuous, in this case we have 

/■0O+27T „ , 

p= / Tr[p$e]$ fl d0. (5-78) 
Using (5-69) Tr[p$ e ] obtain as follows 

Tr[p^ e ] = Tr[p E 6 ™ I ^ >< m \] = 2ir^9 m -!—P PB (9) rn , (5-79) 

m=0 rn S + 1 

where Pp^ is probability of measuring a particular value of phase and is normalized so that 
the integral of Pps^e) over a 2n region of 9 is equal to one. 

P PB ($ e ) = -L V e l(n '-"^ < n | p | n' >=< | p | 9 >, (5-80) 

n,n — 

where thus obtained results are in agreement with those of already obtained by one of the 
authors in [24, 25, 4]. 

A very important subset of these states will be the physical partial phase states, of which 
the coherent state is a particular example. The phase states are themselves unphysical and so 
the best attempt at a physical phase measurement will only project the system into a physical 
partial phase state [24] . In the following, we obtain a physical partial phase state tomography 
i.e., coherent spin states tomography. 
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5.3 Coherent spin states tomography 

To reconstruct a mixed or pure quantum state of a spin s is possible through coherent states: its 
density matrix is fixed by the probabilities to measure the value s along 4s(s+l) appropriately 
chosen directions in space. Thus, after inverting the experimental data, the statistical operator 
is parameterized entirely by expectation values. 

A coherent spin state | n > is associated to each point of the surface of the unit sphere. 

| n >= exp[—i6m((f)).s\ \ s,n z >, (5-81) 

where m(0) = (—sin(p,cos(f),0). 

A stereographic projection of the surface of the sphere to the complex plane give the 

expansion of a coherent state [26] as follows 

1 2s ( \ 1/2 
I s ' n >= 77X1 — ^7 ^ 2s // k z \s-k,n z > . (5-82) 

In order to show that the density matrix p of a spin s is determined unambiguously by ap- 
propriate measurement with a Stern-Gerlach apparatus one precedes as follows. Distribute 
N s = (2s + 1) 2 axes \ s,n > with 1 < n < N s , over (2s+l) cones about the z axis with different 
opening angles such that the set of the (2s+l) directions on each cone is invariant under a 
rotation about z by an angle ^f+i)- 

An unnormalized statistical density operator is then fixed by measuring the N s relative 
frequencies 

Pn{n n ) =< n n | p | n n > , 1 < n < N s , (5-83) 

that is, by the expectation values of the statistical operator p in the coherent states | n n >. 
You obtain N s linear relations between probabilities P n (n n ) and the matrix elements of the 
density matrix with respect to the basis \s — k, n z >. This set of equations can be inverted 
by standard techniques if the directions n n are chosen as described above. For a spin s, the 
projection operators 

I Qn) =1 n n >< n n |, (5-84) 
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constitute thus a quorum Q. In general, a quorum is defined as a collection of (hermitian) 
operators having the property that their expectation values are sufficient to reconstruct the 

quantum state of the system at hand. (Q n | defined as the dual of the quorum (5-84): 

i n s N s 

7^— tt E E I Qn)(Q n ' 1= , 1 < n,n> < N s . (5-85) 

\ ZS l ) n=\n'=\ 

Therefore, this coherent spin state introduced above is same as the phase state. 

In order to obtain spin tomography relation in the finite dimensional Banach space we 
assume that 

p = J Tr(p | Q n ))(Q n | dp n , (5-86) 
is a density matrix in infinite dimensional Banach space. Also let 

p' = E A «l<H (5-87) 

n 

be a density matrix in finite dimensional Banach space which is obtained from truncating the 
infinite dimensional Banach space. Using the properties of density matrix we have 

p - p' > 0, (5-88) 

which is comparison with semidefinite programming and using complementary slackness con- 
dition, we get 

Z(p -p') = or Z(p - \ n | Q n )) = 0. (5-89) 
Similar to supperoperator formalism we obtain 

A n = {Q n | p) = Tr[pQ n ] = P n . (5-90) 

Therefore we obtain the tomography formula in finite dimensional Hilbert space as the follow: 

P S = ^—7H P nQ n , (5-91) 

ZS "T" 1 71=1 



where the coefficients P n satisfy 

< P n < 1 , 1 < n < N s . (5-92) 
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The operators Q n do even define an optimal quorum since exactly (2s + l) 2 numbers have to 
be determined experimentally which equals the number of free real parameters of the (unnor- 
malized) hermitian density matrix p . Thus obtained results are in agreement with those of 
already obtained by one of the authors in [24, 25, 4, 26]. 

It is important to note that, although each of the P n is a probability, they do not sum up 
to unity: 

< ^ P n < (2s + l) 2 (5-93) 

71=1 

. This is due to the fact that they all refer to different orientations of the Stern-Gerlach 
apparatus, being thus associated with the measurement of incompatible observables, 

[Q n ,Qn']^0,l <n,ri < N s , (5-94) 

since the scalar product < n n \n' n > of two coherent states is different from zero. The sum in 
(5-93) cannot take the value (2s + l) 2 since this would require a common eigenstate of all the 
operators Q n which does not exist due to (5-94). By an appropriate choice of the directions 
n n (all in the neighborhood of one single direction no, say), the sum can be arbitrarily close 
to (2s + l) 2 for states peaked about n . Similarly, the sum of all P n cannot take on the value 
zero since this would require a vanishing density matrix which is impossible. If, however, 
considered as a sum of expectation values, there is no need for the numbers P n to sum up 
to unity. Nevertheless, they are not completely independent when arising from a statistical 
operator: its normalization implies that 

i N s 

Tr[p s ] = Tr[— - £ P n Q n ] = 1, (5-95) 

As + 1 n =l 

turning one of the probabilities into a function of the (2s + l) 2 — 1 = 4s(s + 1) others, leaving 
us with the correct number of free real parameters needed to specify a density matrix[26]. 
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Using the elegant method of convex semidefinite optimization method and superoperator for- 
malism, we have been able to obtain the quantum tomography in finite dimensional representa- 
tion for some set of mixed density matrices. In this method we have been able to obtain finite 
qudit, N-qubit quantum tomography, phase tomography and coherent spin state tomography 
, where these results that obtained are in agreement with those of ref[21] and [24, 25, 4, 26] . 
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